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Abstract 


The availability of monoclonal antibodies (mAbs) specific for the SARS-coronavirus (SARS-CoV) is important for the development of 
both diagnostic tools and treatment of infection. A molecular characterization of nine monoclonal antibodies raised in immune mice, using 
highly purified, inactivated SARS-CoV as the inoculating antigen, is presented in this report. These antibodies are specific for numerous viral 
protein targets, and six of them are able to effectively neutralize SARS-CoV in vitro, including one with a neutralizing titre of 0.075 nM. A 
phylogenetic analysis of the heavy and light chain sequences reveals that the mAbs share considerable homology. The majority of the heavy 
chains belong to a single Ig germline V-gene family, while considerably more sequence variation is evident in the light chain sequences. These 
analyses demonstrate that neutralization ability can be correlated with specific murine Vy-gene alleles. For instance, one evident trend is high 
sequence conservation in the Vy chains of the neutralizing mAbs, particularly in CDR-1 and CDR-2. The results suggest that optimization 
of murine mAbs for neutralization of SARS-CoV infection will likely be possible, and will aid in the development of diagnostic tools and 
passive treatments for SARS-CoV infection. 
© 2004 Elsevier Ltd. All rights reserved. 
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1. Introduction infected more than 8000 people in various countries world- 
wide and caused approximately 800 deaths (Drosten et al., 

SARS-coronavirus (SARS-CoV), the causative agent of 2003a,b; WHO; _ http://www.who.int/csr/sars/country/en/). 
severe acute respiratory syndrome (SARS) in humans, has The whole genomes of SARS-CoV isolates, implicated in 
the 2003 outbreak in Toronto, have been sequenced and char- 
acterized (Marra et al., 2003; Rota et al., 2003). Characteri- 
zation of this virus continues at a phenomenal rate, and our 
understanding of the function of numerous viral proteins, the 


Abbreviations: CDR, complementarity determining region; CoV, coro- 
navirus; D, diversity; F.C.A., Freund’s complete adjuvant; F.I.A., Freund’s 
incomplete adjuvant; FR, framework; Ig, immunoglobulin; I.P., intraperi- 


toneal; J, joining; PBS, phosphate buffered saline; SARS, severe acute res- phylogeny of SARS-CoV, and the viral life cycle continues 
piratory syndrome; S.C., subcutaneous; Vy, variable region heavy chain; to grow (reviewed in Stadler et al., 2004; Eickmann et al., 
VL, variable region light chain 2003; Thiel et al., 2003). 
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prevention and control of future outbreaks remains the devel- 
opment of passive immunotherapy with SARS-CoV-specific 
antibodies (Holmes, 2003). It has been suggested that pro- 
tection might be afforded by passive immunotherapy with 
concentrated SARS-CoV-specific IgG antibodies (Li et al., 
2003a), and reports have established that infected individ- 
uals can benefit from treatment with serum from recovered 
patients (Pearson et al., 2003). Recently, it was also reported 
that viral replication was inhibited in the lower respiratory 
tract of naive mice if they underwent passive immunization 
with neutralizing antibodies present in immune serum derived 
from infected mice (Subbarao et al., 2004). A monoclonal 
antibody developed from a human non-immune single chain 
variable region (scFv) library has also been shown to neu- 
tralize SARS-CoV in vitro, strengthening the argument that 
passive immunotherapy with highly specific mAbs might be 
very effective in controlling SARS-CoV infection (Sui et al., 
2004). 

The production of mAbs specific for SARS-CoV is vi- 
tal for studies of viral pathogenesis, and the development 
of both diagnostic tools and vaccines. Since the develop- 
ment of serum antibodies after infection with SARS-CoV 
can take 1-3 weeks (Li et al., 2003a), assays that can accu- 
rately detect the presence of viral nucleic acids or proteins 
may be preferred for rapid diagnosis of SARS-CoV infec- 
tion. The profiles of antibody responses to SARS-CoV have 
been well established (Li et al., 2003a). Analyses have also 
identified viral proteins that might serve as the best mark- 
ers for immunological detection of infection by SARS-CoV 
(Tan et al., 2004; Lu et al., 2004). Similarly, the character- 
ization of immunogenic peptides derived from SARS-CoV 
structural proteins has also allowed for the identification of 
epitopes that are recognized by antibodies present in pa- 
tient serum (Wang et al., 2003). With this knowledge, the 
development of mAbs that can be used for both diagnos- 
tic assays and clinical treatments should be an attainable 
goal. 

Herein, we further characterize the immunogenetics and 
neutralizing endpoints of purified murine hybridoma-derived 
mAbs raised in mice, using highly purified SARS-CoV as 
the inoculating antigen. Numerous reports exist, which char- 
acterize antibodies raised against surrogate SARS-CoV im- 
munogens. Examples include polyclonal antibodies raised 
against peptides and a recombinant SARS-CoV nucleopro- 
tein (Chang et al., 2004; Lin et al., 2003), and monoclonal 
antibodies raised against a SARS-CoV-derived recombinant 
protein fragment (Zhou et al., 2004). Similarly, a neutraliz- 
ing mAb specific for the spike protein, derived from naive 
human phage display libraries, has been developed and char- 
acterized (Sui et al., 2004). Little has been published on mAbs 
developed to the native viral particle. The relevant immuno- 
genetic characteristics of a panel of nine murine mAbs raised 
to whole SARS-CoV, six of which can effectively neutralize 
the virus in vitro, are presented. This development provides 
a distinct advantage in the search for an effective passive im- 
mune therapy, because these mAbs were raised against the 


intact virus, rather than any individual viral protein or im- 
munogenic peptide. 


2. Materials and methods 


2.1. Production of highly purified SARS-CoV for 
inoculation of mice 


All of the procedures employed for the production 
of SARS-CoV are discussed in detail in Berry et al. 
(2004). Briefly, SARS-CoV (Tor-3 strain; Krokhin et al., 
2003) was expanded after plaque purification in Vero-6 cell 
monolayers and partially purified through a sucrose cushion, 
and then further purified using iodixanol gradient centrifu- 
gation. Fractions were tested by Western immunoblot with 
convalescent patient serum, and the fractions that reacted 
with SARS-CoV were pooled, dialyzed against phosphate 
buffered saline (PBS), and further concentrated by ultracen- 
trifugation for 1.5h at 150,000 x g, resulting in highly puri- 
fied whole virus particles. 


2.2. Inoculation of mice and production of hybridomas 
and antibodies 


All of the procedures employed for the production of hy- 
bridomas and mAbs are discussed in detail in Berry et al. 
(2004). Briefly, 5-6-week-old female BALB/c mice (Charles 
River, Wilmington, MA) were injected subcutaneous (S.C.) 
with 50g of beta-propiolactone-inactivated SARS-CoV 
(Tor-3 strain) with an equal part of Freund’s complete ad- 
juvant (F.C.A.) (H37-Ra) from Difco (BD, Oakville, ON). 
Thirty days later, the mice received 50 wg of SARS-CoV S.C. 
in Freund’s incomplete adjuvant (F.I.A.) in a total volume of 
100 wl. The mice received 5 wg of the same antigen in a total 
volume of 100 pl S.C. with LA. on days 48 and 63. The mice 
received a final booster inoculation intraperitoneal (I.P.) with 
5 wg of highly purified SARS-CoV in 200 wl PBS 3 days be- 
fore euthanization by anesthesia overdose, spleen removal, 
and subsequent hybridoma fusion. After hybridoma fusion 
and establishment of stable clones, the hybridoma super- 
natants were screened via ELISA using purified SARS-CoV 
as the target antigen, and isotyped with a commercial murine 
isotyping dipstick test (Roche), facilitating the selection of 
appropriate primers for subsequent RT-PCR. Monoclonal an- 
tibodies were produced in medium scale in 450 and 800 cm? 
sterile disposable roller bottle flasks (Corning). Roller bot- 
tles containing 200 ml of complete hybridoma media [BD- 
Quantum yield (VWR, Canada), 10% fetal bovine serum 
(Wisent, Montreal), 10% Biogro-Hybridoma serum-free sup- 
plement (Biogro Technology, Winnipeg)] were inoculated 
with 50 ml of hybridoma cells such that a final concentration 
of more than | million viable cells/ml was achieved. The cul- 
tures were allowed to grow to extinction and were harvested 
at 7-10 days post-inoculation. Supernatants were clarified 
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by removing cell mass by centrifugation, and concentrated 
on Amicon 8400 stirred cell Ultrafiltration-nitrogen concen- 
trators (VWR, Canada) with a YM-30 (Millipore, USA) 
membrane. Individual membranes were used for the con- 
centration of each mAb. The concentrated supernatants were 
mixed 1:1 with protein A-binding buffer (Pierce) and puri- 
fied on a equilibrated 1 or 5 ml protein G-sepharose column 
(Amersham Biotech) as per the manufacturers instructions. 
Eluted antibody was dialyzed to remove salt using Centriprep 
YM-30 Centrifugal Filter Units (Millipore) with a 30kDa 
cutoff. The mAb solutions were filter sterilized through low- 
protein-binding 0.22 wm syringe filters (VWR). Protein con- 
centration was determined using the micro-bca assay (Pierce) 
and the IgG concentration was standardized to 1 mg/ml in 
sterile PBS (Gibco). 


2.3. Virus neutralization assay 


The SARS-CoV in vitro neutralization assay was per- 
formed essentially as described (Berry et al., 2004), with the 
following modifications. The panel of purified mAbs was 
sterile filtered and each sample standardized to a concen- 
tration of | mg/ml of IgG before being subjected to serial 
dilutions and incubation with live SARS-CoV. 


2.4. Cloning of the heavy (Vy) and light (Vi) chains of 
murine mAbs 


Approximately 1-2 million mAb-secreting hybridoma 
cells were collected and homogenized via passage through 
a 20-guage needle using an RNase-free 5 ml syringe. Total 
RNA was isolated using RNeasy® Mini Spin kit according 
to the manufacturer’s instructions (QIAGEN®). cDNA was 
produced via reverse transcription using 10 pg to 10 yg of 
RNA template, a 15 base oligo-dT primer and SuperScript II 
reverse transcriptase (Invitrogen). Reactions were incubated 
at 42 °C for 1h, followed by inactivation of the reverse tran- 
scriptase at 70°C for 15 min. RNA was removed from the 
reactions by incubating them with RNaseH (Invitrogen) at 
37°C for 20 min. The resulting cDNA was used as a tem- 
plate for PCR amplification using Taq polymerase (Invitro- 
gen) and various 5’ and 3’ primers specific for Vy and Vi 
genes and. The primers used are as follows (in 5’—3’ orienta- 
tion), and are described in detail elsewhere (Berry, 1999; Bar- 
bas et al., 2001; Dattamajumdar et al., 1996). Vy-specific, 5’: 
UmIgVH: TGA GGT GCA GCT GGA GGA GTC; MHCL- 
1: ATG GAC TTB GDA YTG AGC T; MHCL-2: ATG GAA 
TGG ASC TGG RTC TTT CTC T; MHCL-3: ATG AAA 
GTG TTG AGT CTG TTG TAC GTG; MHCL-4: ATG RAS 
TTS KGG YTM ARC TKG RTT; MHC9: AGG TI AIC 
TIC TCG AGT CWG G; MSCVH-9: GGT GGT TCC TCT 
AGA TCT TCC CTC GAG GTG AWG YTG GTG GAG 
TC; MSCVH-10: GGT GGT TCC TCT AGA TCT GCC 
CTC GAG GTG CAG SKG GTG GAG TC. Vu-specific, 3’: 
MH125: GGC CAG TGG ATA GAC. V_-specific, 5’: MVL: 
GTG CCA GAT GTG AGC TCG TGA TGA CCC AGT CTC 


CA; MSCVK-9: GGG CCC AGG CGG CCG AGC TCG 
AYA TTG TTC TAC WCC AGT C; MSCVK-10: GGG CCC 
AGG CGG CCG AGC TCG AYA TTG WGC TSA CCC AAT 
C. Vi -specific, 3’: MCK-1: GGA TAC AGT TGG TGC AGC. 
Degenerate nucleotides are indicated as follows: B = T, C, or 
G; D = A, T, or G; I = deoxyinosine; K = Tor G; M=A 
or C; R= AorG;S=CorG; W=AorT; Y=CorT. 
PCR reactions were incubated at 94°C for 2 min, followed 
by 30 cycles of 96°C for 15s, 56°C for 30s, and 72°C for 
2 min. A final extension at 72°C for 10 min completed the 
reactions. Approximately 1—5 wg of the PCR products were 
resolved on 2% agarose gels and stained with ethidium bro- 
mide. PCR products of the appropriate size (350-500 base- 
pairs) were inserted into the Invitrogen pCR2.1 -TOPO-TA™ 
vector and transformed into TOP-10 cells, according to the 
manufacturer’s instructions. Colonies containing the plasmid 
with insert were selected by overnight growth, at 37°C, on 
Luria-Bertani agar plates containing ampicillin (200 pg/ml) 
and X-gal (80 g/ml). White (positive) colonies were incu- 
bated in 5 ml of Luria—Bertani broth containing ampicillin 
(200 g/ml) for 12-16h with vigorous shaking at 37°C, 
and plasmids were isolated using QIAprep® Spin Miniprep 
DNA isolation kits, according to the manufacturer’s instruc- 
tions (QIAGEN®). EcoRI digestion of 1-5 jg of the resulting 
plasmids and resolution of the products via agarose gel elec- 
trophoresis, as described above, confirmed the presence of the 
appropriately sized inserts. One to 10 putative positive clones 
were sequenced with T7 (5’-TAA TAC GAC TCA CTA TAG 
GG-3’) and M13Rev (5’-CAG GAA ACA GCT ATG AC-3’) 
primers, using ABI Prism® BigDye terminator cycle se- 
quencing kits and an ABI Prism® 3100 Genetic Analyzer. 
Consensus sequences were generated using the DNASTAR 
software suite (DNASTAR Inc.). GenBank accession num- 
bers are AY605265 to AY605273 for the Vy sequences, and 
AY605274 to AY605282 for the Vi, sequences presented in 
this work. 


2.5. Identification of Ig germline sequences and 
assignment of relevant regions 


Consensus nucleotide sequences were compared against 
the Mus musculus immunoglobulin (Ig) set database 
using IMGT/V-Quest (Lefranc, 2003; http://imgt.cines. 
fr/home.html). The sequences were concurrently compared 
against the M. musculus Ig germline V-gene database us- 
ing IgBlast (Altschul et al., 1990; http://www.ncbi.nlm.nih. 
gov/igblast/). This allowed the identification of the comple- 
mentarity determining region (CDR) and framework (FR) 
regions of the Vy and Vy sequences, and provided number- 
ing to the inferred amino acid sequences according to Kabat 
et al. (1991). Similarly, the IgBlast results allowed for the 
identification of the most closely related murine Ig germline 
V-genes currently available in these databases. In all cases, 
the entire sequence, including those at the 5’ end of each se- 
quence imposed by the specific primers used in the original 
PCR amplification, were examined. 
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2.6. Sequence alignments 


The inferred amino acid sequences were trimmed to 
remove all residues encoded by the 5’ primer regions, to elim- 
inate potential artefactual matches caused by the primers em- 
ployed in PCR amplification of the heavy and light chains. 
The sequences were then aligned with ClustalW (Thompson 
et al., 1994; http://clustalw.genome.ad.jp/) using standard 
parameters (gap open penalty: 10; gap extension penalty: 
0.05; no weight transition; hydrophilic gaps allowed; weight 
matrix: Blosum). The alignments were then analyzed using 
GeneDoc version 2.6.002 (Nicholas et al., 1997) to produce 
textual alignments and quantify the relatedness of the se- 
quences. 


2.7. Phylogenetic analysis 


Analyses were performed using BioEdit version 5.0.9 
(Hall, 1999) and MEGA version 2.1 software suite (Kumar 
et al., 2001). Only amino acid sequences were examined. 
The sequences, aligned via ClustalW as described above, 
were analyzed via neighbor joining analysis, bootstrapped 
using 1000 replicates. The naive human mAb 80R (Sui et al., 
2004; Vy GenBank accession no. AAS19425; V_ GenBank 
accession no. AAS194372) was initially used as an outlier se- 
quence. This mAb was chosen because it is specific for the 
SARS-CoV spike protein and neutralizes SARS infection in 
vitro, but it was developed from a non-immune human sin- 


Table | 


gle chain variable fragment (scFv) library. A second outlier, 
also human in origin, was found as follows. The Vy, and Vy 
amino acid sequences of F26G18 were compared against the 
human Ig sequences in the NCBI databases using IgBlast, 
as described above. F26G18 was chosen because it was the 
most potent SARS neutralizing murine mAb isolated, and it is 
specific for spike. Human antibodies were chosen to ensure 
that similarity would be minimal, considering that mouse, 
and indeed human, antibodies inherently share considerable 
homology regardless of their specificity. The sequences that 
showed the lowest homology to the Vy and Vy, sequences 
of F26G18 were chosen as follows: VL: GenBank acces- 
sion no. BACO1727, level of identity: 83/115 amino acids; 
Vu: GenBank accession no. CAD19025, level of identity: 
83/132 amino acids. In all of the phylogenetic analyses, the 
human sequences clearly served as effective outliers, and only 
BAC01727 and CAD19025 are included in the relevant fig- 
ures. 


3. Results 


3.1. Summary of the relevant genetic properties of the 
mAbs specific for SARS-CoV 


Table 1 summarizes the relevant characteristics of the 
murine mAbs generated against highly purified, inactivated 
SARS-CoV. Examining the Vy chains, only two Ig V-gene 


Summary of relevant properties of the Vy and Vy; regions of murine mAbs specific for SARS-CoV 


mAb* Class GenBank Target” Neutralizing J-gene® D-gene* Most closely related Ig Percent identity to most closely 
accession titre (nM) germline V-gene4 related Ig germline V-gene° 
numbers 
Nucleotide Amino acid 
Vu 
F26G1 G2a/k AY605271 Spike - Ju4 SP2.6/7/8 VOx-1 92 85 
F26G6 G2b/k = AY605272 Spike - Ju2 None J558.33 95 89 
F26G8 G2a/k AY605273 Spike - Ju2 None J558.33 95 86 
F26G18 G2b/k =AY605269 Spike 0.075 Ju4 FL16.2 J558.50 97 91 
F26G19 G2a/k AY605270 Spike 1 Ju4 Q52.01 J558.50 97 89 
F26G3 G2a/k AY605265 U 26 Jul FL16.1 J558.5 96 94 
F26G7 G2b/k = AY605266 U 6 Ju4 Q52.01 J558.50 97 91 
F26G9 G2a/k AY605267 U 1 Jul SP2.3/4/5 J558.50 97 94 
F26G10 G2a/k AY605268 U 1 Jul SP2.3/4/5 J558.50 96 91 
VL 
F26G1 G2a/k AY605280 Spike - Jq2 N.A. ay4 99 96 
F26G6 G2b/k  =AY605281 Spike - Jx4 N.A. bw20 94 86 
F26G8 G2a/k AY605282 Spike - Jx4 N.A. bw20 92 81 
F26G18 G2b/k = AY605278 Spike 0.075 Jq2 N.A. ce9 97 92 
F26G19 G2a/k AY605279 Spike 1 Jl N.A. cw9 98 96 
F26G3 G2a/k AY605274 U 26 Jq2 N.A. 21-7 98 92 
F26G7 G2b/k = AY605275 U 6 Jl N.A. cw9 98 95 
F26G9 G2a/k AY605276 U 1 Jq2 N.A. 19-25 98 94 
F26G10 G2a/k AY605277 U 1 Jq2 N.A. 19-25 99 98 


* Neutralizing mAbs are indicated by bold text. 
> U: unknown. 


© Most closely related genes, as determined by V-quest and NCBI IgBlast analysis of nucleotide sequences. N.A.: not applicable. 


4 As determined by NCBI IgBlast analysis of nucleotide sequences. 


° Expressed as percentage of identical residues (pairwise), as determined by NCBI IgBlast of nucleotide or inferred amino acid sequences, as indicated. 
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(A) 


J558.50: TGGGACTGTGCTGGCAGGGCCTGGGGCTTCAGTGAAGATGTCCTGCAAGGCTTCTGGCTACACCTTTACCAGCTACTGGATGCACTGGGTAAAACAGAGGCCTGGACAGGGTCTGGAATGGATTGGC 


F26G18: 
F26G19: 
F26G7 : 
F26G9 : 
F26G10: 


J558.50: 
F26G18: 
F26G19: 
F26G7 : 
F26G9 : 
F26G10: 


J558.50: 
F26G18: 
F26G19: 
F26G7 : 
F26G9 : 
F26G10: 


(B) FRI CDR-1 FR2 CDR-2 


J558.50: SGTVLAGPGASVKMSCKASGYTFT SYWMH WVKQRPGQGLEWIG ATIYPGNSDTSYNQKFKG KAKLTAVTSTSTAYMELSSLTNEDSAVYYCTR 
s 


F26G18: ..... PRE stars) akeheaieierace een a e-S- aude 


F26G19: wane IRB seer iia Bie Gi armaresaaraaye ESRe.< 
E26G/ of aguas Rewseveees Be cseuee Tewee 
F26G9 2 eee Reweescerccenes See teen 
F26G1LO: ..4..5 Risie iid wees tee ates eu S23 wasae 


bub eherepaes Sts Nesetersceness, <i ReMlast a age: devatie 


Fig. 1. Alignments of the nucleotide and amino acid sequences of Vy regions of the neutralizing mAbs belonging to the J558 V-gene family. (A) Nucleotide 
alignment of the Vy regions of the mAbs most closely related to the J558.50 V-gene. The consensus sequence used is that of J558.50 (nucleotide sequence 
GenBank accession no. AF303881). A dot in the individual sequences denotes nucleotides that are the same as the consensus. A dash in the individual sequences 
denotes a deletion. Neutralizing clones are shown in bold text. (B) Amino acid alignment of the Vy regions of the mAbs most closely related to the J558.50 Ig 
V-gene. The consensus sequence used is that of J558.50 (amino acid sequence GenBank accession no. AAG39162). A dot in the individual sequences denotes 
amino acids that are the same as the consensus. A dash in the individual sequences denotes a deletion. The framework and complementarity determining regions 


(CDR) are indicated above the appropriate sequence segments in the figure. 


families are represented. Eight out of the nine mAbs share 
high sequence identity with the J558 V-gene family, with the 
remaining mAb, F26G1, sharing high sequence identity with 
the VOx-1 gene. Of the J558 V-gene family members, five 
of the eight mAbs are most closely related to the J558.50 
V-gene, and this set of five mAbs neutralizes SARS-CoV in 
vitro. mAbs F26G6, F26G8 and F26G3 are most closely re- 
lated to the J558.33 and J558.5 V-genes, respectively, which 
are themselves closely related to J558.50. A closer analysis of 
the Vy chains of mAbs sharing significant homology with the 
J558.50 V-gene is presented in Fig. | and Table 2. An align- 
ment of the nucleotide sequences of these Vy sequences re- 
veals a very high level of identity (Fig. 1A). Similarly, at the 
amino acid level, these sequences share considerable iden- 
tity both amongst themselves and with the J558.50 V-gene 
(Fig. 1B), with some mutation away from the germline se- 
quence, which is typical of immunoglobulin V-gene cDNAs 


Table 2 
Comparison of the Vy regions of the SARS-CoV-specific neutralizing mAbs 
sharing significant homology with the J558.50 Ig germline V-gene 


mAb Level of identity vs. J558.50 

Nucleotide Amino acid 
F26G18 267/274 (97%) 86/92 (93%) 
F26G19 266/274 (97%) 84/92 (91%) 
F26G7 267/274 (97%) 87/92 (95%) 
F26G9 267/274 (97%) 89/92 (97%) 
F26G10 265/274 (97%) 86/92 (93%) 


Level of identity (pairwise) is expressed as number of identical residues/total 
number of residues and percent identical residues. 


from T-dependant B cell responses (Fish and Manser, 1987). 
A quantification of the level of identity of these sequences 
is summarized in Table 2. An analysis of the J-genes shows 
that three Jy genes are represented in the Vy chains, at a 
relatively even distribution. There appears to be a slight bias 
towards Jy4 (four mAbs) with Jy 1 (three mAbs) and Jy2 (two 
mAbs) also represented (Table 1). Analysis of the CDR-3 se- 
quences reveals that Dy-regions are present in seven out of 
nine of the mAbs, with homology to numerous Dy-gene se- 
quences available in the databases. A detailed analysis of the 
Dy-regions is presented in the following section. 

Examination of the Vj, chains demonstrates that a wider 
variety of V-genes is represented. Six distinct V-genes are 
present, which suggests that a larger pool of V-genes can be 
selected, compared to the Vy chains, to produce mAbs that 
are specific for SARS-CoV. As in the case for Vy, three J- 
gene families are represented, with a definite bias towards 
Jx2 (five mAbs) and an even distribution of J, 1 and J,4 (two 
mAbs each). 


3.2. Analysis of the Vy and Vi, sequences of the mAb 
panel 


An alignment of all of the amino acid sequences of the Vy 
and Vy, chains is presented in Fig. 2A. An inspection of the 
Vu sequences reveals that CDR-1 and CDR-?2 share a rela- 
tively high level of identity in all of the mAbs. Conversely, 
CDR-3 shows a considerably higher level of sequence vari- 
ation. A close inspection of the Vy sequences shows that 
F26G6/F26G8, F26G7/F26G19, and F26G9/F26G10 are all 
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Fig. 2. (A) Amino acid alignment of the Vy and Vj regions of the murine mAbs raised against SARS-CoV. The consensus sequence (top line) shows amino acid 
residues that are conserved in >50% of the sequences. A dot in the consensus sequence denotes amino acids that are conserved in <50% of the listed sequences. 
A dot in the individual sequences denotes amino acids that are the same as the consensus. A dash in the individual sequences denotes a deletion. Neutralizing 
mAbs are shown in bold text. The framework (FR) and complementarity determining regions (CDR) are indicated above the appropriate sequences in the 
figure. (B) Nucleotide alignments of identifiable D-regions in the Vy CDR3 regions. The most closely related D-regions were determined by V-quest analysis 
as outlined in Section 2. Dots represent identical residues, while dashes represent gaps or nucleotides that were not taken into account for the alignments. 


Neutralizing mAbs are shown in bold text. 


closely related (within each pair), with very few amino acid 
changes evident between these sequences. 

An examination of the Vy sequences reveals that unlike 
the case of Vy, CDR-1, -2, and -3 exhibit considerably higher 
levels of sequence variation. As is the case for the Vy chains, 
the Vy; chains of the F26G6/F26G8, F26G7/F26G19, and 
F26G9/F26G10 pairs share a high level of sequence identity 
(within each pair) and are therefore likely clonally related. 
In the case of both the Vy and Vy, sequences, mAb F26G1 
stands out as being the least closely related to any of the other 
sequences in this panel of mAbs. A quantitative analysis of 
sequence relatedness is discussed in the subsequent sections. 

The discernible Dy-regions present in the CDR-3 of the 
Vu domains show considerable variation, and alignments 
of these regions with available related Dy-region sequences 
identified by V-quest analysis are shown in Fig. 2B. mAb 
F26G1 contains the most well delineated but shortest Dy- 
region, with high homology to SP2.6/7/8 Dy-genes over six 
nucleotides. mAb F26G3 shows high homology to FL16.1 
and considerable homology with SP2.7/8/9 Dy-genes, which 
also overlaps with FL16.1. mAbs F26G7 and F26G19 con- 
tain identical Dy-regions, which share homology with both 
SP2.5/7/8 and Q52.01/02, although delineation of the re- 
gions of homology is less clear. Similarly, mAbs F26G9 
and F26G10 contain identical Dy-regions, with homology 
to SP2.3/4/5, DP2.1, and FL16.2. Manual alignment of the 


sequences suggests that the central portion of the Dy region 
of F26G18 is most likely encoded by an SP2 or DP2 ele- 
ment, with considerable junctional diversity (D) on either end 
(Fig. 2B). mAb F26G18 contains a Dy-region with homol- 
ogy to SP2.5/7, DP2.1, and FL16.2. Discernable Dy-regions 
are lacking in mAbs F26G6 and F26G8, and these mAbs 
are therefore not included in this figure. It is clear from the 
analysis that the Dy-regions of all of these mAbs are highly 
variable, and likely resulted from complex somatic gene re- 
arrangements and junctional diversification during assembly. 


3.3. Quantitative analysis of the relatedness of the 
regions encoded within the Vy and Vi, chains 


Examining the full length of the Vy sequences, as outlined 
in Fig. 2A, a clear pattern emerges regarding amino acid se- 
quence identity. The most closely related Vy sequences are 
F26G9/F26G10 (97% identical), F26G7/F26G19 (95% iden- 
tical), F26G6/F26G8 (92% identical) and F26G18/F26G19 
(82% identical). Both Vy CDR-1 and CDR-2 exhibit high 
levels of identity, with all of the neutralizing mAbs showing 
remarkable sequence conservation, while Vy CDR-3 shows 
a high level of sequence diversity. One neutralizing mAb 
F26G3 shares little identity with any of the other neutralizing 
mAbs because it is derived from a different germline Vy al- 
lele (see the following section for a more detailed analysis 
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Fig. 3. Summary of the percentage identity of the SARS-CoV-specific 
mAbs. (A) Summary of the percentage of pairwise identity, at the amino 
acid level (Y-axis), between the Vy chains of various categories of mAbs 
(all, neutralizing, non-neutralizing) as indicated on the X-axis. White bars: 
full length of sequences; hatched bars: CDR-1; grey bars: CDR-2; black 
bars: CDR-3. (B) Summary of the percentage of identity at the amino acid 
level of the Vi sequences. This figure is labeled the same as (A). 


of this mAb). The most closely related Vj, sequences are 
F26G7/F26G 19 (98% identical), F26G9/F26G10 (96% iden- 
tical), and F26G6/F26G8 (90% identical). The V_. CDRs ex- 
hibit more sequence variability than those found in the Vy 
chains, and as was the case for F26G3 Vy, the Vy, chain 
of this mAb shows very little sequence identity with the Vy 
chains of the other neutralizing mAbs (Fig. 2A). For both 
the Vy and Vy chains, mAb F26G1 demonstrates the lowest 
sequence identity compared to any of the other mAbs in the 
panel (Fig. 2A and B). 

The average level of identity of the heavy and light chains, 
subdivided into pertinent regions, is shown in Fig. 3. Across 
the full length of Vu, it is evident that the average level of 
identity is highest amongst the neutralizing mAbs, albeit by 
a small margin (Fig. 3A). An examination of the Vy CDRs 
reveals that CDR-1 and CDR-2 exhibit the highest level of 
sequence identity, whether subdivided into groups contain- 
ing all mAbs, neutralizing mAbs, or non-neutralizing mAbs. 
However, this finding is most apparent when examining the 
neutralizing mAbs, whose CDR-1 and CDR-2 regions share 
73% and 82% average sequence identity, respectively. On 
the other hand, CDR-3 exhibits the highest sequence vari- 
ability, whether grouped into all mAbs, neutralizing mAbs, 
or non-neutralizing mAbs, with an average level of identity 
ranging from 24% to 47%. An extended CDR-3 domain is 
characteristic of the Vy chains of the SARS-CoV neutraliz- 
ing mAbs. This domain of the neutralizing mAbs contains, 
on an average, 12 amino acid residues, while the average 
amongst all of the mAbs is 10 residues (not including D- 


element negative clones F26G6 and G8). The CDR-3 regions 
of the non-neutralizing mAbs contain only five residues. The 
V_ regions exhibit a lower average level of sequence identity, 
as outlined in Fig. 3B. As with the Vy sequences, the highest 
average level of identity across the full length of the Vi se- 
quences is exhibited by the neutralizing mAbs. However, the 
average level of sequence identity is lower overall in the Vy 
sequences than in the corresponding Vy sequences (Fig. 2A 
and compare Fig. 3A and B). Unlike the Vy sequences, there 
is no clear bias towards a higher level of identity in the 
CDRs, with all three CDRs sharing relatively low average se- 
quence identity, regardless of the grouping of the mAbs into 
all, neutralizing, or non-neutralizing categories. This finding 
again suggests that the V_ CDRs likely influence the neu- 
tralization ability of the mAbs to a lesser extent than the Vy 
CDRs. 


3.4. Summary of the phylogenetic relationship of the 
mAbs 


The overall relatedness of the Vy and Vi sequences of the 
mAbs can be quickly summarized via phylogenetic neighbor 
joining analysis. Fig. 4 shows phylogenetic trees generated 
by an examination of the amino acid sequences of the Vy and 
Vi chains. Fig. 4A summarizes the Vy chains. Clearly, the 
neutralizing mAbs F26G7, F26G9, F26G10, F26G18, and 
F26G19 cluster together and are closely related. Amongst 
these Vy chains, F26G18 is the least closely related to the 
remaining members of this group. Interestingly, F26G18 and 
F26G19, which are both spike-specific neutralizing mAbs, 
cluster further apart, and are in fact on separate branches of 
the tree. Another neutralizing mAb, F26G3, clusters closer to 
mAbs F26G6 and F26G8, which are non-neutralizing, spike- 
specific mAbs. Unfortunately, the protein target of mAb 
F26G3 is currently unknown. However, one would not be 
surprised if this mAb does in fact target spike, but this inter- 
action has been undetectable by the methods employed thus 
far (Berry et al., 2004). 

Although these analyses cannot determine whether 
F26G6/F26G8, F26G7/F26G19, and F26G9/F26G10 are in 
fact mAbs derived from the same clones, it is clear that they 
cluster in related branches according to their overall proper- 
ties (e.g. neutralizing, spike-specific). Itis also clear that mAb 
F26G1 is least related to any of the other mAbs in this panel, 
occupying its own branch well separated from the rest of the 
mAbs. Examining the Vy chains sharing the highest iden- 
tity with the J558.50 V-gene, all of which are neutralizing, 
it is evident that a similar pattern emerges. Monoclonal anti- 
bodies F26G9/F26G10 and F26G7/F26G19 are encoded by 
V-genes that cluster closely together, with F26G18 occupying 
a divergent branch on the tree (Fig. 4B). Essentially the same 
pattern emerges when examining the Vy, chains (Fig. 4C), 
with neutralizing and spike-specific mAbs clustering in re- 
lated branches, and F26G1 occupying a separate branch on 
the tree, illustrating that this mAb is the least related to any 
of the other mAbs in the panel. 
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Fig. 4. Phylogenetic analysis of the mAbs specific for SARS-CoV. Trees are based on neighbor joining analysis of the inferred amino acid sequences of the 
relevant Vy and Vi, sequences, as outlined in Section 2. (A) Vy regions of all mAbs. (B) Vy regions of all mAbs sharing significant identity with the J558.50 
Ig germline V-gene. (C) Vy, regions of all mAbs. In each case, the outliner sequences are of human origin, chosen as outlined in Section 2. Neutralizing mAbs 
are shown in bold text. The number at each node represents the level of bootstrap support (expressed as a percentage) for the node over the total number of 


replicates performed. Only those values above 50% are reported. 


3.5. Detailed analysis of the neutralizing mAbs 


Table 3 presents a summary of the level of identity evi- 
dent in the group of mAbs that exhibit SARS-CoV neutraliz- 
ing ability, relative to the most efficiently neutralizing mAb, 
F26G18. As shown in this table, the neutralizing titre varies 
from a low of 0.075 nM (F26G18) to a high almost 350-fold 
greater, at 26 nM (F26G3). As a comparison, it was recently 
reported that a human mAb, 80R, developed from a naive 
human immune scFv library and specific for the angiotensin- 


converting enzyme 2 (ACE2) binding domain of spike, ex- 
hibited a neutralizing concentration as low as 0.37 nM (Sui 
et al., 2004). The relationship between F26G18 and the rest 
of the neutralizing mAbs was therefore examined as an at- 
tempt to identify any common theme or relationship that ex- 
ists amongst them. Clearly, mAb F26G3 is least related to 
F26G18 or any of the other neutralizing mAbs, in both the 
Vy and Vy chains (Table 3). Examining the Vy chains, the 
average level of identity (compared to F26G18) for all of the 
neutralizing mAbs is high across the entire sequence of the 


Table 3 

Comparison of the Vy and Vj regions of the most efficiently neutralizing mAb (F26G18) with the corresponding regions of the remainder of the neutralizing 

mAbs 

mAb Neutralizing titre (nM) Region examined and level of identity vs. F26G18* 

Full length CDR-1 CDR-2 CDR-3 

Vu 
F26G19 1 95/116 (82%) 3/5 (60%) 15/17 (88%) 3/13 (23%) 
F26G9 1 96/116 (83%) 5/5 (100%) 16/17 (94%) 2/13 (15%) 
F26G10 1 93/116 (80%) 5/5 (100%) 15/17 (88%) 2/13 (15%) 
F26G7 6 99/116 (93%) 4/5 (80%) 16/17 (94%) 4/13 (31%) 
F26G3 26 77/116 (66%) 3/5 (60%) 11/17 (64%) 4/13 (31%) 
Average identity’ (%) 81+ 10 80 + 20 86 + 12 2348 

VL 
F26G19 1 69/99 (70%) 8/11 (72%) 2/7 (28%) 3/9 (33%) 
F26G9 1 65/99 (66%) 5/11 (45%) 2/7 (28%) 7/9 (17%) 
F26G10 1 64/99 (65%) 5/11 (63%) 2/7 (28%) 6/9 (66%) 
F26G7 6 70/99 (71%) 8/11 (72%) 2/7 (28%) 3/9 (33%) 
F26G3 26 59/103 (57%) S/15 (33%) 3/7 (42%) 2/9 (22%) 
Average identity’ (%) 66 + 6 57+17 3146 46 + 24 


* Level of identity (pairwise) is expressed as number of identical residues/total number of residues and percent identical residues. F26G18 is the most 


efficiently neutralizing mAb, with a neutralizing titre of 0.075 nM. 


> The average level of identity (+S.D.) is calculated from the values of the specified regions for each given mAb. 
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Vy chain, and within CDR-1 and CDR-2. Conversely, CDR- 
3 is the least closely related region in all of the neutralizing 
mAbs. The same pattern is evident when the average level of 
sequence identity of the Vy chains across all mAb sequences 
is examined (Fig. 3A). Examining the Vy, chains, the average 
level of identity between the rest of the neutralizing mAbs 
and F26G18 is relatively low (Table 3), which is similar to 
the relatively low average level of identity exhibited by all 
of the mAbs in the panel (Fig. 3B). The average level of 
identity in CDR-1 and CDR-2 of Vy. is lower, compared to 
F26G18, than the corresponding regions in Vy, while CDR-3 
exhibits a higher average level of identity. In the Vy, chains, 
a slightly lower average level of identity is evident in CDR- 
2, compared to CDR-1 and CDR-3; however, this is far less 
noticeable than the difference in relatedness between CDR-3 
compared to CDR-1 and CDR-2 in the Vy chains. 


4. Discussion and conclusions 


This panel of mAbs derived from mice immunized with 
highly purified whole SARS-CoV exhibits diverse specificity 
for binding targets and variable in vitro neutralization ability. 
Neutralization titres amongst the mAbs vary from the low- 
est, at 0.075 nM for F26G18, to a high of 26 nM, for F26G3. 
Interestingly, two of the mAbs that exhibit the highest neu- 
tralizing efficiency, F26G18 and F26G19, are specific for the 
spike protein (Berry et al., 2004). As has been established 
conclusively, spike is required for binding of SARS-CoV to 
its receptor, ACE2, and mAbs that block this interaction can 
neutralize SARS-CoV (Li et al., 2003b; Sui et al., 2004). 
These results suggest that F26G18 and F26G19 likely neu- 
tralize SARS-CoV infection by directly targeting spike, al- 
though whether the interaction between spike and its receptor 
is targeted specifically remains to be determined, and ongo- 
ing work to answer this question is in progress. Although 
five of the mAbs are confirmed to be spike-specific, three of 
these (F26G1, F26G6, and F26G8) are unable to neutralize 
SARS-CoV infection in vitro. These results suggest that this 
panel of mAbs recognizes a diverse set of epitopes within the 
spike protein, some of which are in regions of spike that are 
not involved in viral adhesion or other steps in the infection 
process. Of the panel of six neutralizing mAbs, the targets 
of four have not yet been identified, and work is currently 
underway to determine them. 

Examining the level of pairwise identity of the neutraliz- 
ing mAbs at the amino acid level reveals that the highest se- 
quence variability is exhibited by the Vy CDR-3, while, not 
surprisingly, CDR-1 and CDR-2 of the Vy chains demon- 
strate a greater level of sequence identity. Conversely, all 
three CDRs of the Vy chains exhibit similar levels of se- 
quence identity, which are, on average, lower than the levels 
of identity in the corresponding Vy regions. The variability 
in CDR-3 of the Vy chains is therefore likely a significant 
contributing factor to the target specificity of these mAbs, 
and the Vi CDRs likely influence the neutralization abil- 


ity of the mAbs to a lesser extent than the Vy CDRs. The 
average length of CDR-3 also appears to influence neutraliz- 
ing ability. The Vy CDR-3 of the neutralizing mAbs contains 
more amino acids than the average number in the whole panel 
of mAbs, while the CDR-3 regions of the non-neutralizing 
mAbs contain fewer residues. Conversely, the average length 
of Vi CDR-3 is constant. These results suggest that both the 
length and sequence of the Vy CDR-3 influence the neutral- 
ization ability of the mAbs, most likely by determining the 
specificity of contact residues for key neutralizing epitopes 
on SARS-CoV proteins. 

Within the group of neutralizing mAbs, a trend is evi- 
dent that relates homology in the Vy chains of the most po- 
tent neutralizing mAb, F26G18, with those that neutralize 
SARS-CoV less efficiently. The clustering of the sequences of 
the neutralizing mAbs based on their relatedness to F26G18 
tends to correlate well with their ability to neutralize SARS- 
CoV in vitro. For both Vy and VL, the mAb with the lowest 
neutralizing efficiency, F26G3, has the lowest overall se- 
quence identity when compared to mAb F26G18, and when 
compared to the average sequence identity across the whole 
panel of neutralizing mAbs. This trend is more apparent for 
the Vy sequences, since the percentage of identity of the 
VL chains relative to F26G18 is lower overall. This observa- 
tion is consistent with the dominant role of Vy compared to 
V{ in previously identified neutralizing mAbs against a viral 
pathogen (Barbas et al., 1993). CDR-1 and CDR-?2 of the Vy 
chains of the neutralizing mAbs are most closely related to 
the same regions in F26G18, while CDR-3 is less conserved. 
From this analysis, one can infer that high homology in Vy 
CDR-1 and CDR-2 (relative to F26G18) of the neutralizing 
mAbs appears to correlate with neutralizing ability. Inter- 
estingly, F26G7 exhibits a six-fold higher neutralizing titre 
than F26G19, even though their sequences are closely related 
both to each other and to F26G18, revealing that the trend for 
relatedness to F26G18 correlating with neutralizing ability 
does not follow for all of the mAbs in this panel. This ob- 
servation also suggests that small changes in the amino acid 
composition of these mAbs can significantly alter neutral- 
ization efficiency. Overall, this analysis clearly demonstrates 
that the Vy regions of the neutralizing mAbs are more highly 
related than the Vy; chains, and suggests that the sequence of 
the Vi, chains plays a smaller role in neutralizing ability than 
the Vy chains. 

Analysis of the sequences of the mAbs reveals that the 
Vu and Vy chains of several mAbs are closely related. The 
pairs F26G6/F26G8, F26G7/F26G19, and F26G9/G10 ex- 
hibit few differences at the amino acid level, within each 
pair. This observation could be explained in several ways. 
First, each pair might have originated from the same B cell 
progenitor and they are clonally related (Fish et al., 1989), 
with the few observed sequence variations resulting from 
somatic mutation. Second, identical CDR-3 rearrangements 
might have occurred in the B cell progenitors, and antigen se- 
lection may have produced several independently rearranged 
but closely related B cell clones. Third and alternatively, 
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sequence changes might have occurred during reverse tran- 
scription or PCR amplification, and the mAbs originated 
from the same, independently picked hybridoma clone. While 
these analyses do not allow a differentiation between these 
potential causes of this high sequence similarity, we would 
speculate that the first explanation is most likely true. 

The Vy chains of eight of the mAbs belong to the J558 
V-gene family, and six of these mAbs neutralize SARS-CoV 
in vitro. This result suggests that a small pool of B cells, 
expressing related Ig V-genes, was selected for in response 
to exposure to SARS-CoV. This is not entirely surprising, 
considering that the J558 V-gene family is estimated to con- 
tain several hundred members in the genome of the BALB/c 
mouse (Livant et al., 1986) and this family is highly selected 
in murine B cells (Haines et al., 2001). It is clear that a re- 
stricted pool of Vy genes was selected to encode SARS-CoV 
neutralizing mAbs. A bias towards the selection of genes in 
the J558 V-gene family to generate autoantibodies in the tight- 
skin mouse has also been reported (Kasturi et al., 1994). Jy 
gene usage is distributed relatively evenly in the mAbs gen- 
erated during this work, which agrees well with the observa- 
tion that Jy gene usage in several murine innate antibodies 
belonging to the J558 V-gene family, is approximately evenly 
distributed (Seidl et al., 1999). These results suggest that a 
typical variety of joining (J) genes is selected during the gen- 
eration of Vy chains in these clones. Conversely, a bias was 
evident in the selection of J genes in the generation of the 
V, chains in these clones, which agrees with published re- 
ports of preferential selection of J, 1 and J,2 genes in adult 
mouse splenic cells (Wood and Coleclough, 1984; Nishi et 
al., 1985). This observation suggests that the Vy, chains of 
these mAbs appear to be derived from a more diverse pool 
of V genes, and may be less critical for contributing contact 
residues, relative to Vy, for binding to SARS-CoV proteins. 
A complete absence of Dy-regions, and as a result a small 
CDR-3 in the V-domain, in two of the non-neutralizing, spike- 
specific mAbs, F26G6 and F26G8, suggests that this region 
contributes significantly to neutralization. The diversity of 
the length, sequences, and apparent order of assembly of the 
Dy-regions present in these mAbs suggests that significant 
somatic alteration of sequences flanking the Dy-regions oc- 
curred during junctional assembly, or that somatic mutations 
occurred during affinity maturation in a germinal center. The 
presence of identical Dy elements in the CDR-3 regions of 
F26G7/F26G19 and F26G9/F26G 10 supports the notion that 
these pairs are clonally related, and derived from the same 
parental V-gene rearrangement. 

In summary, multiple mAbs able to neutralize SARS-CoV 
infection in vitro were generated in this work, and they dis- 
play neutralization titres in the nanomolar to picomolar range. 
The majority of the Vy chains of the neutralizing mAbs share 
significant homology to the J558 Ig V-gene family, suggesting 
selection of a limited pool of V-genes occurred in mice ex- 
posed to highly purified SARS-CoV. High levels of sequence 
identity in CDR-1 and CDR-?2 of the Vy chains, compared to 
low levels of sequence identity in all of the Vy chains, sug- 


gests that neutralization ability is imparted to these mAbs by 
contact residues within the Vy chains. Chain shuffling exper- 
iments will be employed to determine the relative contribu- 
tions of the Vy and Vy chains to SARS-CoV neutralization 
(Berry et al., 2003). 

Passive immune therapy is a viable treatment for expo- 
sure to numerous viruses. The development of the mAbs 
reported here is an important first step in isolating mAbs gen- 
erated against the intact SARS virus that could be used for 
passive immunotherapy. Humanization of murine mAbs is a 
viable method for the development of potential passive im- 
munotherapeutics (Roguska et al., 1994; Rader et al., 1998), 
and these methods can be applied to the mAbs developed in 
this work. One caveat to this analysis is that it is not an ex- 
haustive statistical analysis of the relationships between the 
mAbs described here. Due to the small numbers of mAbs 
examined, only trends and general relationships regarding 
the levels of sequence identity between any of the mAbs and 
their immunological properties can be inferred. Regardless, 
this analysis provides a useful summary of murine mAbs that 
possess the potential for development into viable immune 
therapies and diagnostic tools for SARS-CoV infection. The 
recurrent usage of particular murine V-gene elements to en- 
code virus-neutralizing antibodies is becoming more evident 
and detailed molecular analyses become available (Kalinke 
et al., 1996). Studies with animal models are currently un- 
derway to determine the in vivo neutralization properties of 
the neutralizing mAbs summarized here. 
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